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Abstract Chloroplasts regulate their growth to optimize photosynthesis. Quantitative data 
shows that the ratio of total chloroplast area to mesophyll cell area is constant across dif¬ 
ferent cells within a single species, and also across species. Wild-type chloroplasts exhibit 
little scatter around this trend; highly irregularly-shaped mutant chloroplasts exhibit more 
scatter. Here we propose a model motivated by a bacterial quorum-sensing model consist¬ 
ing of a switch-like signalling network that turns off chloroplast growth. We calculated the 
dependence of the location of the relevant saddle-node bifurcation on the geometry of the 
chloroplasts. Our model exhibits a linear trend, with linearly growing scatter dependent on 
chloroplast shape, consistent with the data. When modelled chloroplasts are of a shape that 
grows with a constant area to volume ratio (disks, cylinders) we find a linear trend with min¬ 
imal scatter. Chloroplasts with area and volume that do not grow proportionally (spheres) 
exhibit a linear trend with additional scatter. 

Keywords chloroplasts • ordinary differential equations • bifurcation ■ model ■ growth 
regulation • switch 


1 Introduction 


Chloroplasts are the organelles within plant cells which convert sunlight to usable energy 
via photosynthesis. Their exposure to sunlight is integral to optimal photosynthetic output, 
so their coverage of a cell’s surface area is important. There is evidence that chloroplast 
growth and replication is regulated to control the percentage of a mesophyll cell’s surface 
area that is occupied by chloroplasts. This is quantified by comparing the plan areas of the 
chloroplasts and the cell, defined as the areas taken up by the chloroplasts and the cell when 
viewed from above ( |Pyke and Leech[[l991| l. Several studies support this claim. 

[Ellis and Leech| ( l985| l found that leaves of Triticum aestivum L. and T. monocoecum L. 
demonstrated a trade-off between chloroplast size and number that resulted in a similar total 
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Fig. 1 Figure reprinted (with permission) from 
ship between total chloroplast plan area per cell and cell plan area for each of the wild type and three mutants 
of Arabidopsis thaliana. The slope of all four lines is similar. There are two layers of chloroplasts within a 
cell, one at the top and one at the bottom, so “full” cells would exhibit a slope of 2. Scatter varies, with the 
wild-type and arcl mutants having the least scatter {R^ — 0.92 and — 0.93, respectively, in (a) and (b)) 
and arcl and arc3 mutants showing more scatter {R^ — 0.81 and ^ — 0.61, respectively, in (c) and (d)). 


chloroplast plan area per unit cell plan area. |Pyke and Leech| \\992) then examined meso- 
phyll cells from a wild type and three chloroplast replication and growth mutants (termed 
“arc” mutants) of Arabidopsis thaliana and consistent with the observations of Ellis and 
Leech, they found an inverse relationship between the mean individual chloroplast area and 
the total number of chloroplasts (see Figure [TJ. For each of the wild-type and the three mu¬ 
tants, the relationship between total chloroplast plan area and cell plan area was linear, no 
matter the individual size or shape of the chloroplasts. Moreover, the slope of the data plotted 
in the “plan area plane” (total chloroplast plan area versus cell plan area) was quantitatively 
similar for all four data sets. Two of the mutants showed considerably more scatter than the 
wild type and the other mutant, and the mutant with the most scatter was observed to have 
irregularly shaped chloroplasts. 

|Osteryoung et al| l |l998| t unearthed more evidence of this phenomenon from a collection 
of Arabidopsis cells of both the wild type and two transgenic types. Their figure plotting the 
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figure 5, showing the chloroplast 


Fig. 2 Figure reprinted (with permission) from 
plan areas versus cell size for wild type (WT, open circles) and two different transgenic types (1-1, filled 
boxes, and 2-1, filled triangles) of Arabidopsis. All three relationships are linear with similar slopes; see 
calculated slopes and values. 


total chloroplast plan area versus the cell plan area for each cell is duplicated here in Figure 
Similar to the data in |Pyke and Leech| ^1992\ , each of the three types displayed linear 
relationships of similar slope, with scatter that scaled with cell size. The transgenic types 
- with larger and more irregularly shaped chloroplasts than the wild type - presented with 
more scatter; see the values in Figure]^ 

|Pyke| ( |1999| l further showed that this regulated chloroplast density persists not only 
across different cells, but across different species. He found a linear relationship between 
total chloroplast area and cell plan area when considering an average data point from each 
of eight different species and several Arabidopsis mutants. He noted that the signalling meth¬ 
ods used by chloroplasts to regulate their growth and division are not known, but concluded 
that “chloroplast density in relation to the size of the cell seems to be involved”. 

The most notable feature in the data sets shown in Figures [T] and is the linear rela¬ 
tionship between total chloroplast plan area and cell plan area; this is backed up by the data 
from |Pyk^ ( |1999| l. The second notable feature is that the scatter about the relationship grows 
somewhat linearly with cell plan area, more like a percentage of the plan area ratio rather 
than an absolute measure. This data can be resummarized in the context of the ratio of plan 
areas, which we’ll call a, defined as the total chloroplast plan area divided by the cell plan 
area. These two observations - the linearity of the relationship and the linearly increasing 
scatter - can be reduced down to the observation that the range of ratios of plan areas is 
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constant for all cell plan areas. The third important observation is that irregularly shaped 
chloroplasts have more scatter, which means that chloroplast geometry matters. 

Here we propose a mechanism to explain this observed chloroplast growth regulation 
phenomenon. The ability to sense the collective density of one’s own population is found 
in certain bacteria and is referred to as “quorum sensing” l |Dockery and Keener] [2001^. 
Much work has been focused on understanding quorum sensing; Dockery and Keen^ 2001^ 
modelled the biochemical quorum-sensing switch in Pseudomonas Aeruginosa. 

At the core of our model is the Dockery-Keener biochemical switch. In bacteria, flip¬ 
ping this switch triggers a shift in gene regulation allowing a sufficiently large population to 
collectively change behaviour. These triggered behaviours include producing exopolysac¬ 
charide to form a biofilm, secreting toxins and turning on bioluminescence. In chloroplasts, 
we propose that a similar switch turns off chloroplast growth. 

Although chloroplasts may have a quorum-sensing biochemical pathway homologous 
to that of bacteria, we prefer to think of this quorum-sensing model simply as a biochem¬ 
ical system with switch-like behaviour. This switch-like behaviour is the important feature 
for our results. Even if the chloroplast growth regulation pathway operates differently, any 
switch-like behaviour involving secretion of the signal ought to behave similarly to what we 
describe here. We find that the model, when applied to a system of chloroplasts within a cell, 
can capture both the linear trend and the linearly growing scatter in the data. In the model, 
chloroplast shape influences the extent of the scatter. 


2 Mathematical model 

2.1 Cell and chloroplast geometry 

Quorum sensing is dependent on the geometry of the system. In a system comprised of 
chloroplasts within a leaf cell, the biochemical processes of secretion, degradation, and ab¬ 
sorption are coupled to geometry via dependence on the surface area of the chloroplasts and 
the volume density of the chloroplasts within the cell. We define the geometry of the sys¬ 
tem in order to find explicit definitions of the chloroplast density, p, and the chloroplasts’ 
effective membrane permeability, 5 - which is dependent on surface area - to use in the 
quorum-sensing model. 

Consider a cell of fixed size that contains a population of chloroplasts. A leaf mesophyll 
cell of A. thaliana is largely occupied by one or more vacuoles, leaving only a thin layer of 
the cytosol available for the chloroplasts at the surface ( |Pyke[ [2009111999| ; s^e Figure[^ We 
model this available cytosolic space as two thin slices that can each hold a single layer of 
chloroplasts, with a top surface area of S and a thickness of T. The volume of a cell is thus 
V = 2St, and its plan area (area viewed from above) is S. We look at a range of cell sizes by 
considering a range of S values; we keep the thickness of each cytosolic space constant at T. 

Chloroplast shape is more complicated. Chloroplasts typically resemble flattened foot¬ 
balls, although this varies widely l |Pyke[[2009| l. The variation in their morphology is both 
characteristic of different plant species or cell types, and dynamic in response to the environ¬ 
ment. For example, the arc6 mutant of A. thaliana has highly irregularly shaped chloroplasts 
in comparison to the wild-type chloroplasts ( jPyke and Leech[[l994| l. Dynamically, chloro¬ 
plasts can deform in response to being squeezed against neighbouring plastids or vacuoles 
or the cell wall, which can happen over long time scales as the cell grows, or in shorter 
time scales as chloroplasts move in response to changes in light ( |Pyke| [2009] [Briggs and] 
[Christie! [2002[l. Chloroplast morphology also changes during division when the dividing 
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Fig. 3 This cartoon summarizes the features of A. thaliana leaf cell images (see examples in Pyke et al 
(l99^ and [Hall and Langdale|jT99^ ) from both (a) an above view and (b) a transverse section view. The 
large vacuole occupying the middle of the cell pushes the chloroplasts into a thin layer against the cell wall, 
usually concentrated on the top and bottom for better access to light sources. Based on this cartoon, we model 
the relevant cytosolic space as two thin layers of thickness T and top surface area S that can each hold a single 
layer of chloroplasts, one at the top and one at the bottom. Figure created in OmniGraffie. 



Fig. 4 Chloroplasts are modelled either as all cylinders as in (a); all disks as in (b); or all spheres as in (c). 
Orientation is as shown, within a layer of fixed thickness T. Each shape grows in the direction labelled f; 
all other dimensions are fixed. Our model contains two such layers, of top surface area S. Figure created in 
Adobe Illustrator. 


chloroplast goes through a “dumbell-shaped” stage l [Pyke[[2009| t. Since there is such vari¬ 
ability in chloroplast shape, we represent chloroplasts as either all cylinders, all thin disks, 
or all spheres, and we fix all the dimensions except for one which we refer to as the growth 
direction, 1. Cylinders are of a fixed radius t/2 but variable length disks are a fixed thick¬ 
ness T with variable radius t, spheres have a variable radius I < t/2. See Figure]^ 

For any of these three basic chloroplast shapes, we can define the surface area of each 
chloroplast as Sc, and the volume of each chloroplast as Vc- In a cell of volume V = 2St, the 
chloroplast density is then p = nVc(2ST and the effective permeability is 5 = nScB, where 
n is the number of chloroplasts and 5 is the per surface area permeability of the chloroplast 
membrane. 
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2.2 Quorum sensing equations 

In this section, we reformalize the Dockery-Keener model to make it suitable for the chloro- 
plast context by explicitly incorporating the parameters that govern cell and chloroplast 
geometry (p,5). 

[Dockery and Keen^ ( |200l] l based their quorum sensing model on the eight component 
gene-regulatory system described in [Van Delden and Iglewski] ( |1998[ l. They simplified it 
down to two ordinary differential equations using time scale assumptions. We adopt the two 
variable model in this paper but consider chloroplasts within a cell in place of cells within a 
matrix, and chloroplast density instead of cell density. The first equation tracks the concen¬ 
tration of an autoinducer (A), a signalling molecule which is secreted and can be sensed by 
all chloroplasts. The second equation tracks the concentration of a protein (f?) that dimerizes 
with the autoinducer to form an activator of autoinducer production, thereby generating a 
positive feedback loop. When the autoinducer reaches a critical level, we assume (but do not 
explicitly model) that a downstream pathway is triggered and chloroplast growth is shut off. 
The equations for A and R are 


where 


dR RA „ , „ 

(1) 

dA RA 

dt Ka+RA^^'^ 

(2) 

-p) 

p kE{\-p) + 8 

(3) 


Here kE and are degradation rates of the signalling molecule outside and inside the 
chloroplast, respectively; kn is the degradation rate of the protein R inside the chloroplast; 
Vr, Va control the maximum rate of Michaelis-Menten-type production, and Ka, Kr control 
the shape of the production curve. Aq and Rq are background production rates. The parameter 
d depends on p, the relative volume density of chloroplasts (with maximum 1), and 5, the 
effective permeability of the chloroplast membrane. This complicated expression arises from 
the reduction of the system from eight variables down to two. It consists of two terms. The 
first term is the degradation rate constant for autoinducer inside the chloroplast. The second 
term gives the secretion rate when the cytosolic concentration of autoinducer is assumed to 
be in quasi-steady state. We use d as the parameter for our bifurcation analysis because it is 
the aggregate parameter through which the system depends on chloroplast geometry (5) and 
population density (p). Using the geometry defined in Section [ZT| the bifurcation parameter 
d can be written as a function of f and n for each fixed cell surface S (denoted by a subscript 

S): 


ds{i,n) =kA + 


2SrSc{t)5 

kE{l-"^)+nSS)S’ 


(4) 


where Sc and W will vary depending on the shape and size of the model chloroplasts. When 
d is written in this form it becomes clear that chloroplast geometry and growth plays an 
important role in the dynamics of this system. 

The steady states are found by solving a cubic in RA. Figure]^ shows a bifurcation plot 
of the steady state value of A versus d. Each of the two “knees” indicates a saddle-node 
bifurcation. As the chloroplasts grow, d decreases, so the biologically relevant bifurcation is 
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Fig. 5 The solid curve is a bifurcation plot showing the concentration of signalling molecule A versus the bi¬ 
furcation parameter d. A cell with low chloroplast density will begin on the right side of the lower branch and 
move left (d will decrease) as the chloroplasts grow. When d reaches the bifurcation point aXd*, the system 
will fall off the knee and the concentration of signalling molecule A will rise above the critical value A* (indi¬ 
cated by a horizontal dashed line). This transition is the end-of-growth bifurcation at which the chloroplasts 
will stop growing. The bifurcation diagram was calculated using AUTO and plotted using MATLAB. 


the one at the lower d value; we label this J-value as J*. A cell with low chloroplast density 
p will begin on the lower branch at the right-hand side (the low steady state) and then track 
left as the chloroplasts grow, falling off the knee onto the upper branch (the high steady 
state). This bifurcation - depicted in the phase plane in Figure]^- is a biochemical switch. 

We assume that the chloroplasts’ growth is binary (on or off) dependent on the magni¬ 
tude of A. While A is below some critical value A*, they remain in the growth regime. When 
A rises above A*, the chloroplasts stop growing. We refer to this switch as the end-of-growth 
bifurcation. 

We are interested in determining the relationship between size, number, and shape of 
chloroplasts at the critical d* value, and how well that relationship predicts the data pre¬ 
sented by |Pyke and Leech| l |1992| l and [Osteryoung et al| ( |1998| l. In particular we aim to 
determine how the linear trend and linearly growing scatter of the data, and the variance in 
the amount of scatter, depend on chloroplast geometry. 


3 Analysis 

3.1 General procedure 

The following analysis is performed for each of the three chloroplast geometries. 

We fix cell size S and denote dependence on cell size by a subscript S where relevant. 
Since the chloroplasts vary in only one dimension the bifurcation parameter ds is a function 
of only two variables: the growth dimension, and the number of chloroplasts n. The form 
varies depending on which geometry we consider. Our goal is to find which combinations 
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Fig. 6 The nullclines in the phase plane of the signalling molecule A vs the protein R, plotted for three values 
of the bifurcation parameter d, show the passage through the first saddle-node bifurcation (the end-of-growth 
bifurcation). As chloroplast density increases in a cell of fixed size, d will begin above d* and decrease 
through the bifurcation point at d*. The curve A' = 0 is solid; S' = 0 is dashed, and is invariant to changes in 
d. Figure created in MATLAB. 


of size and chloroplast number correspond to the bifurcation point at ds = d*, and then 
calculate the corresponding plan areas of those size-number combinations. Using those data 
points we can then draw the plan area plane (total chloroplast plan area versus cell plan area 
S) for each case. 

Finding the bifurcation point requires solving the equation ds{£,n) = d*. We nondimen- 
sionalize this equation and reformulate it as a search for the zeros of the function gs{£,n) 
defined by 


gs{e,n) = {D* - K) 


SMt: 



2St ) 


nVdt)\ I nSc(t)' 
2Sz J S 


(5) 


where D* is the nondimensional bifurcation point defined as D* = d*/kE, Ps = ^e/SS, and 
K = k^/kE- Each (£,n) pair corresponding to the bifurcation point describes a size £ and 
number n of chloroplasts at which growth shuts off. The allowable pairs are restricted by 
the physical requirement that the chloroplasts must fit within the cell: ”'^ 2 Sz ^ Note that 
n is discrete, so the bifurcation “curve” in the £ — n plane consists of a collection of discrete 
points. 

Once we have a collection of these (f, n) pairs that correspond to the bifurcation curve, 
we can calculate the ratio of the total chloroplast plan area to the cell plan area - the slope 
of the line through the origin and data point in the plan area plane - for each pair. We label 
this ratio generally as as{£,n) \ an explicit definition depends on the chloroplast geometry 
that we’re working with. We then plot the points (5, tty • S) in the plan area plane. 

We repeat these calculations over a range of fixed S values to generate the full data set. 

The scatter and linearity of the data can also be examined analytically for each geometric 
case by looking at the expressions for cty; we include this below on a case by case basis. The 
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range of as values dictates the vertical scatter of the data at that fixed S'-value on the x-axis. 
It is possible that different {i,n) pairs return the same as value so the vertical scatter does 
not necessarily correlate to the number of bifurcation points. The variance in the range of as 
values across the different cell sizes S dictates the shape or trend of the data as a whole. If 
the range of values of as is invariant in S (i.e., as = a), then the predicted chloroplast plan 
area to cell plan area relationship is linear with linearly growing scatter. 


3.2 Cylinder-shaped chloroplasts 


A right cylinder of variable length £ and fixed diameter T has a surface area of Sc = jl -£ 
nx£ and volume Vc = To keep the chloroplast dimensions in a physically reasonable 

regime we consider i> t/2. When the chloroplasts take this shape in a cell of fixed surface 
area S and thickness T, the function in ([^ becomes 


gs{l,n) = {D* 



(1- 

85 / 


n7txi \ 

85 ) 


j_ Ttttn / T j_ 1 \ 
5 V 2/ ^ / 


( 6 ) 


To find the roots - the {£, n) pairs that correspond to the end-of-growth bifurcation for this S 
- we set gs{£, n) = 0. This can be solved numerically using Newton’s method (for the general 
case), or in this case analytically as a quadratic equation in i (taking only the positive root). 
Then the plan area ratio for cylinders, defined as 

x£n 

^^cylinder ^ 


is calculated for each {£, n) pair at the end-of-growth bifurcation. Each unique value 

corresponds to a point {S, a • S') in the plan area plane. The range of a ■ S values provides the 
vertical distribution of the data at that cell size S. 

We repeat this process for a range of S values, generating a set of data points at each S. A 
sample plot showing the scatter of plan areas for the full distribution of predicted chloroplast 
size-number pairs for cylinder-shaped chloroplasts at the end-of-growth bifurcation is shown 
in Figure]^ The apparent slope of the model data (S, a ■ S) depends on as yet unconstrained 
parameter values. In the figures, we have chosen kA so as to get a slope close to what is 
observed in mesophyll cells. All parameter choices are listed in Appendix [A| 

To gain further insight into the scatter and linear trend of this distribution, we solve 0 
for n and replace n in = 0 from Q, rearranging to find an expression for in terms 

^ 8 8(l + ^)-fo(D*-j^) 

Scylind.r ;j. g ( J ^ ^ ^ ^ ^ 1 

This is subject to the restrictions D* — K > 0 and jis > 0, and we require a > 0. To examine 
the scatter of points at each S value, consider varying i. Under the domain constraint l> xjT., 
this expression has little variability as a function of I and thus only a small range of a values 
for each S. This translates to the small amount of scatter that we see in Figure Next, to 
examine the overall trend of the data as cell size S varies, we note that the only ^-dependence 
is in fis = kE/S5. Since S is large with respect to Re/ 5, jis is small for any S. This means 
that there is little variance in with respect to S, which explains the linear trend seen 

in Figure I?] To see both of these observations more clearly, we use the fact that fis is small 
to express O^Scyimd„ leading order in jis as 


8 


( 9 ) 
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Fig. 7 A sample of calculated cylindrical chloroplast total plan areas that correspond to the end-of-growth 
bifurcation, over a range of cell plan areas. Note the linear trend, with little scatter. Areas are in Figure 
created in MATLAB. [Parameters: S = 5, ks = 0.1, = O.4.] 


This leading order is constant in both i and S, which predicts a linear trend with 

linearly growing but minimal scatter. These analytical observations match what we observe 
in our calculated data in Figure|7] 


3.3 Disk-shaped and spherical chloroplasts 

We perform the above analysis for disk-shaped chloroplasts, with a fixed thickness of T and 
a variable radius of i, and spherical chloroplasts with a variable radius i restricted to be 
£ < t/2. The disk shape has a surface area of Sc = 27t£{£ -F t) and a volume Vc = the 
sphere has surface area Sc = 47ri'^ and volume Vc = Both disks and spheres look the 
same from above, resulting in identical plan area ratio definitions: 

Tzfin Tzfin 

~ ^ ’ ^S^,pkere ~ ^ 

Again, we explicitly calculate the (^,n) pairs that correspond to the bifurcation point by 
finding the roots of ([^ with the appropriate Sdi) and V’e(£) expressions, and then calculate 
the as values for each (£,n) pair. We repeat this for each cell size S, plotting all {S, as ■ S) 
points in the plan area plane. See this calculated data for disk-shaped chloroplasts in Figure 
[^and spherical chloroplasts in Figure]^ 

To explain the observed scatter and the linear trend in these figures we look at the ex¬ 
pressions for as- As with the cylinders, we use the definition of as in each case ( |10[ l and 
each geometry’s version of gj = 0 from l|^ to write as as a function of only i, analogous to 
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The calculated expression for a as a function of I for the disk is 


4(l + |)-fo(D*-j^) 


( 11 ) 


The disk expression has slightly more variance with t than the cylinder due to the factor of 
(1 + 1) instead of the cylinder’s (1 + ^); this shows up as slightly more scatter in the data. 
The scatter is still minimal, though, which is evidenced by the leading order expression 


2 


( 12 ) 


which is independent of i. As in the cylinder case, the scatter in the disk case is also restricted 
to the higher order terms. The sphere expression, on the other hand, is 

etc = 3 t 6^-MD*-K) 
s.pher, 2e6j{D*-K+l)-l3s{D*-K) ^ 

which is similar in form to the other two but multiplied by |. That factor appears due to the 
way the 5 and p factors simplify. It translates to a larger range of as values for each cell size 
S due to the stronger dependence on i, which means that the spherical chloroplasts exhibit 
considerably more scatter than the cylindrical or disk-shaped chloroplasts. 

The leading-order expression for the sphere as is 


3t 1 

2£ (D*_^:+l)’ 


(14) 


from which we clearly see that there is f!-dependent scatter even at leading order. However, 
®+p;isre independent of S at leading order, which agrees with the linear trend of the 

data. 


4 Discussion 

Here we’ve proposed a model to understand the linear trend and scatter seen in the data from 
|Pyke and Leech| ( [1992[ l and |Osteryoung et~n| ( |1998l l. Our model explains both the linearity 
and scatter as being dependent on the geometry of growth. For cylinders or disks that grow 
in one dimension, our model predicts a linear trend with minimal scatter. For spheres, we see 
more scatter. This is consistent with mutant phenotypes that are characterized by irregular 
shapes and greater scatter when compared to the more cylindrical wild-type chloroplasts. 

4.1 Geometry effects 

If the geometry of a chloroplast has a surface area, plan area, and volume that scale together 
as it grows, at the end-of-growth bifurcation we get a narrow range of plan area ratios a for 
each cell size for any reasonable range of parameter values. For example, the cylindrical- 
and disk-shaped chloroplasts, whose surface area, plan area, and volume all scale with £ and 
£^, respectively, show minimal scatter in the model. If the areas and volume do not scale 
together, then we get a wide range of plan area ratios for each cell size, which shows up as 
scatter. This appears in the case of spherical chloroplasts, which have a volume that scales 
with £^ but surface and plan areas that scale with In all cases, whether the range of a 
values is small or large, the actual a values have little dependence on cell size S. As a result 
we observe a linear relationship between total chloroplast plan area and cell plan area, with 
linearly growing scatter. 
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Fig. 8 A sample of calculated disk-shaped chloroplast total plan areas that correspond to the end-of-growth 
bifurcation, over a range of cell plan areas. Note the linear trend, with some scatter. Areas are in Figure 
created in MATLAB. [Parameters: 5 = 5, = 0.1, = 0.45.] 
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Fig. 9 A sample of calculated spherical chloropla.st total plan areas that correspond to the end-of-growth 
bifurcation, over a range of cell plan areas. Note the linear trend, with significant scatter. Areas are in flm^. 
Figure created in MATLAB. [Parameters: 5 = 5, = 0.1, kn = 0.25.] 







A model of chloroplast growth regulation in mesophyll cells 


13 


4.2 Data comparison 


The values for the data presented by Pyke and Leech (1992 1 indicate that although all 
the data sets have a linear trend, the amount of scatter varies amongst the wild-type and 
the mutants (see Figure [T}. Our model can explain the variance in scatter by considering 
geometry. The data sets with little scatter - the wild type, with = 0.92, and the arcl 
mutant with R^ = 0.93 - can be explained by chloroplasts that have areas and volume that 
scale together as they grow, such as the cylinders or disks examined here. The data sets with 
additional scatter - the arcl mutant, with R?' = 0.81, and the arc3 mutant with R^ = 0.61 - 
can be explained by chloroplasts that have areas and volumes that scale disproportionally, 
such as spheres. All three shapes considered here produce linear trends as seen in the data. 

Revisiting the original data and measuring either the number of chloroplasts in each cell 
or the size of the chloroplasts, for each data point, would allow us to predict where those 
points should fall in the plan area plane based on our model. A comparison of the residuals 
for a least squares fit could tell us whether our model accounts better for the scatter in the 
data than noise would. 

The features of the data from |Osteryoung et al| (|1998|l are comparable to those from| Pyk^ 
[and Leech| ( |1992| l, and can be similarly explained by our model if the transgenic types with 
higher scatter have chloroplast shapes with volume and area that scale disproportionally. 


4.3 Model feasibility 

Our model proposes that chloroplasts regulate their growth via a geometry-dependent bio¬ 
chemical switch. This idea came directly out of the bacterial quorum-sensing literature with¬ 
out justification. We do note that the endosymbionic theory traces the origins of chloroplasts 
and mitochondria to bacteria ( [Martin and Herrmann[ |1998| l. Moreover, jPykej ( |1999| ( com¬ 
ments that “the discovery of bacterial cell division gene homologs strongly suggests that 
higher plant plastid division is based on a system that has evolved from that utilized in 
prokaryotic cells,” leaving open the possibility that plastid (i.e., chloroplast) density regula¬ 
tion could also have evolved from a bacterial quorum-sensing system. 

That possibility aside, our model suggests constraints on what features are required of 
a system similar to what we’ve proposed: a biochemically-triggered growth switch, and 
some chloroplast-surface-area-dependent secreted signal which is produced proportional to 
chloroplast volume. This latter requirement could be complicated depending on where the 
genes reside (nuclear or chloroplastular). 


A Parameter ehoiees 

The ODE system constants in equations ^ and were taken directly from [Dockery and KeenerH2001k 
Vs = 2, Va = 2, A'r = 1, kt = 1, Ro = 0.05, Aq = 0.05, = 0.7. With these constants, the relevant bifurcation 

value of d is at d* = 0.5568. Units of time and concentration have been implicitly nondimensionalized. 

The parameters of cell thickness T and surface area S, rates and and permeability 5, that determine 

the value of the bifurcation parameter d are chosen as follows. _ 

Cell size S is variable but of order 10^ (see Fi gure [T] as duplicated from jPyke and Leech|jl99^ ) 
and cell thickness T is constant and of order 1 /im jPyke|[T999k We take T = 4 pm. 

The expression for d requires d* —k^ > 0; this provides an acceptable range for choosing kjt, < 0.55. 
Since k/^ and represent the degradation of the autoinducer (inside and outside the chloroplasts, respec¬ 
tively), we assume they are of the same order of magnitude. In addition, we need to choose 5 such that 
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S5 » kE\ this is the regime in which the autoinducer can collect to a sufficient density outside the chloro- 
plasts to trigger the switch in d. Final specification of kE and 5 is subject to a rough fit of the observed 
slope of the data in |Pyke and Leech] j 1992) to the leading order expression for a for each geometry. We 
choose to fix kE — 0.1 and 5 = 5 for all geometries, and alter k^ to adjust the slope for each geometry. There 
is some flexibility in choosing these parameters within an order of magnitude or more. 
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